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RESEARCH  MEMORANDUM 

PREDICTION  OF  FLAME  VELOCITIES  OF  HYDROCARBON  FLAMES 
By  Gordon  L-  Dugger  and  Dorothy  M.  Simon 


SUMMARY 

The  laminar -flame -velocity  datu  previously  reported  by  the  Lewis 
laboratory  are  surveyed  with  respect  to  the  correspondence  between 
experimental  flame  velocities  and  values  predicted  by  semltheoreticai 
and  empirical  methods-  The  combustible  mixture  variables  covered  are 
hydrocarbon  structure  (56  hydrocarbons),  equivalence  ratio  of  fuel-air 
mixture,  mole  fraction  of  oxygen  in  the  primary  oxygen-nitrogen  mixture 
(0-17  to  0-50),  and  initial  mixture  temperature  (200°  to  615°  K) • The 
semltheoreticai  methods  of  prediction  considered  are  based  on  three 
approximate  theoretical  equations  for  flame  velocity:  the  Semenov 

equation  (.bimolecular) , which  is  primarily  based  on  conductive  heat 
transfer  between  the  flame  and  the  reactants;  the  Tanford -Pease  equa- 
tion, which  is  based  on  the  diffusion  of  chain  carriers  of  the  oxidation 
reaction  into  the  reactants;  and  the  Manson  equation,  which  is  a modi- 
fication of  the  monentum-pressure-drop  equation  that  does  not  include 
chemical  kinetics-  In  euch  equation  a semiempirical  factor  is  used  to 
bring  the  predicted  values  for  a given  variable  and  fuel  into  the  best 
average  agreement  with  the  data,  so  that  the  variation  in  the  relative 
prediction  of  the  individual  datum  points  may  be  considered.  For  the 
resulting  semltheoreticai  equations,  it  is  assumed  that  thermal  equilib- 
rium is  attained  at,  the  end  of  the  flame  zone,  and  vulues  for  the  trans- 
port properties  are  estimated  by  extrapolation  and  simple  additive  rela- 
tions. The  empirical  relations  between  combustible  mixture  variables 
and  flame  velocity  are  based  on  the  usual  methods  of  correlation. 

Within  these  limitations,  the  results  may  be  generalized  as  follows: 

1.  The  three  semltheoreticai  equations  predict  relative  flame 
velocities  reasonably  well,  generally  with  mean  deviations  of  2 to 
15  percent. 

2.  Considering  only  the  bimolecular  fuel-oxygen  reaction  and  with 
low-temperature  activation  energies,  the  Semenov  equation  may  be  used  to 
give  good  relative  predictions  when  a "sterlc  factor"  is  determined 
semiempirically . 

3.  The  average  "specific  rate  constants"  (kg)  obtained  with  the 
Tanford-Pease  equation  for  molecular  fuel  - active  particle  reactions 
behave  as  Arrhenius  rate  constants  only  with  the  data  on  the  effect  of 
initial  temperature,  not  with  the  equivalence-ratio  or  oxygen- 
concentration  data.  The  relative  predictions  obtained  by  this  equation 
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are  not  very  sensitive  either  to  the  temperature  dependence  assigned  to 
the  diffusion  coefficients  or  to  the  recombination  factor  computed  for 
hydrogen  atoms.  Regarding  active  particles  considered,  better  results 
are  generally  obtained  for  hydrocarbon  flames  when  the  hydroxyl  radical 
and  the  hydrogen  and  oxygen  atoms  (OH,  H,  and  0)  are  considered  rather 
than  H alone;  this  is  particularly  true  for  studies  of  the  effect  of 
equivalence  ratio- 

4-  For  the  equivalence- ratio  and  oxygen-concentration  data,  con- 
siderably better  relative  predictions  are  obtained  from  the  Manson 
equation  when  the  pressure  drop  across  the  flame  front  is  considered  to 
be  due  to  H,  OH,  and  0 than  when  only  H atoms  are  taken  into  account- 

5-  The  example  presented  by  the  ethylene  data  shows  that  an  empiri- 
cally determined  rate  constant  or  proportionality  factor  from  equivalence - 
ratio  data  at  atmospheric  pressure  and  room  temperature  may  predict  the 
effect  of  initial  temperature  or  oxygen  concentration  within  approxi- 
mately 20  percent  by  the  Semenov  or  Tanford-Pease  equation  and  within 
approximately  30  percent  by  the  Manson  equation - 

6-  For" engineering  applications,  the  effects  of  the  parameters 
studied  could  be  estimated  just  as  satisfactorily,  and  more  easily,  by 
one  or  another  of  the  ec^irical  correlations  indicated,  as  compared  with 
the  three  semi theoretical  equations  considered-  However,  the  use  of  the 
semitheoretical  equations  in  some  cases  reduces  the  number  of  constants 
required. 


INTRODUCTION 

The  ability  to  predict  flame  velocities  of  fuels  is  of  growing 
importance  in  the  field  of  aircraft  propulsion,  since  a correlation  has 
been  found  between  combustion  efficiency  of  a ram-jet  burner  and  the 
laminar  flame  velocity  of  the  fuel  (reference  1).  The  prediction  of 
flame  velocities  is  difficult  for  three  reasons,  (l)  There  is  no  com- 
plete, rigorous  theory  which  can  be  readily  applied-  There  are,  however, 
a number  of  approximate  equations  in  the  literature  which  approach  the 
problem  of  flame  propagation  from  various  viewpoints.  (2)  There  are  no 
data  on  the  kinetics  of  the  oxidation  process  under  flame  conditions 
and  very  few  data  on  transport  properties  at  high  temperatures.  (3) 

Different  methods  of  flame -velocity  measurement  give  different  values, 
so  that  it  is  difficult  to  compare  data  from  different  sources-  The 
uncertainty  in  measurements  made  by  a given  method  is  of  the  order  of 
5 percent. 

* 

In  this  report  flame-velocity  measurements  made  at  the  NACA  (ref- 
erences 2 to  7)  for  different  hydrocarbons,  initial  temperatures,  and 
compositions  are  used  with  semitheoretical  and  empirical  methods  of  t 
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flame-velocity  prediction  to  show  the  correspondence  between  the  meas- 
ured velocities  and  the  predicted  velocities.  The  semitheoretical 
methods  are  based  on  the  Semenov  equation  (reference  8),  the  Tanford- 
Pease  squai-e-root  law  (reference  9),  and  the  Manson  equation  (refer- 
ence ll) - These  three  equations  were  derived  from  different  models  of 
the  flame -propagation  process-  The  Semenov  model  is  essentially  a 
thermal  model  which  includes  chemical  reaction  kinetics;  the  Tanford- 
Pease  model  is  based  on  the  diffusion  of  chain  carriers  of  the  oxidation 
reaction;  and  Manson  used  a modification  of  the  momentum-pressure-drop 
equation  which  does  not  Include  chemical  kinetics-  The  empirical  rela- 
tions  for  the  effect  on  maximum  flame  speed  of  hydrocarbon  structure,  of 
initial  mixture  temperature,  and  of  oxygen  concentration  were  based  on 
the  usual  methods  of  correlation- 

It  is  recognized  that  there  are  other  approximate  theoretical  equa- 
tions which  may  give  as  good  or  better  predictions  of  the  NACA  experi- 
mental results.  These  particular  equations  were  chosen  because  they 
exemplify  three  different  approaches  to  the  problem  and  because  the 
Semenov  and  Tanford-Pease  equations  were  used  in  previous  NACA  papers 
and  many  of  the  calculations  had  already  been  made  individually.  It  is 
also  recognized  that  there  are  many  flame -speed  data  in  the  literature 
whicn  could  be  used  in  such  studies;  the  present  paper  is  confined  to 
NACA  data  because  many  of  the  calculations  had  been  made,  and  because  it 
was  desirable  to  avoid  the  complications  arising  from  comparing  data 
obtained  by  many  different  methods  and  techniques. 


NOMENCLATURE 

The  following  nomenclature  is  used  in  this  paper: 
a fuel  concentration,  molecules  of  fuel  per  cm3  of  mixture 

b oxygen  concentration,  molecules  of  oxygen  per  cm3  of  mixture 

B-j_  term  near  unity  arising  frcsn  radical  recombination 
C0  molar  heat  capacity  at  constant  pressure,  cal/mole-°K 

Cp  specific  heat,  cal/g-°K 

Cp  mean  specific  heat,  Tq  to  Tf,  cal/g-°K 

D diffusion  coefficient,  cm2/sec 

Dj_  diffusion  coefficient  of  ith  active  species  into  unburned  gas  at 
initial  mixture  temperature,  cm2/ sec 
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E activation  energy,  kcal/g-mole 

fa  fitted  activation  energy  for  a group  of  data  points,  kcal/g-mcle 

K empirical  proportionality  constant  between  experimental  flame 

velocity  and  value  predicted  by  the  Manson  equation 

T0 

Kg  empirical  K from  Manson  equation  when  Ap  *»  l/2  pg 

empirical  K from  Manson  equation  when 

a ^lo(  D°h  D0\T0 

Ap  = 1/2  ^pH  + PQg  + % Dg/  Tf 

Kg,Kj-  average  Kg  or  K£  for  a group  of  data  points 

k specific  rate  constant,  cm^/molecule-sec 

kg  weighted  mean  kj  for  the  three  active  particles  H,  OH,  and 

0,  each  reacting  with  fuel  molecules,  cm^/molecule-sec 

1^  average  for  a group  of  data  points,  cm^/molecule-sec 

kp  specific  rate  constant  for  reaction  between  fuel  molecule  and 

pth  active  particle,  cur/molecule-sec 

Ljjj  total  concentration  of  gas  at  mean  combustion- zone  temperature, 

molecule/cnr^ 

M molecular  weight 

m molecularity  of  flame  reaction 

mole  fraction  of  1th  component 

n total  number  of  molecules  of  HgO  and  CC>2  in  products  per  molecule 

of  fuel  by  stoichiometric  relation 

n^/ rv,  moles  of  reactants  per  mole  of  products  from  stoichiometric 

equation 

P steric  factor,  or  probability  factor,  from  the  expression 

k - 

P average  P for  a group  of  data  points 

p total  pressure  of  mixture 
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Pj  mole  fraction  of  the  active  particle  in  burned  gas 
Q mole  fraction  of  potential  combustion  product  in  unburned  gas 

Q'  mole  fraction  of  fuel  in  unburned  gas 

R universal  gas  constant,  kcal/(g-mole) (°K)  or  ergs/(g-mole) (°K) 

T absolute  temperature,  °K 

computed  equilibrium  flame  temperature,  °K 

U flame  velocity,  cm/sec;  implies  maximum  flame  velocity  with 
respect  to  equivalence  ratio,  except  in  discussion  of  the 
variation  of  equivalence  ratio 

v reaction  rate,  molecules  reacting  per  cm^/sec 

Z collision  number;  the  number  of  molecular  collisions  per  second 
when  the  concentration  is  one  molecule  of  each  type  per  cm3 

a mole  fraction  of  oxygen  in  oxygen-nitrogen  portion  of  mixture, 
°2/^°2+N2^  ' 

6 percentage  mean  deviation  in  the  ratio  of  predicted  flame  velocity 
to  the  experimental  value  for  a given  group  of  data 

q viscosity  of  mixture,  poise 

tj ^ viscosity  of  itb  component,  poise 

6m  ratio  of  mean  reaction- zone  temperature  to  initial  temperature 
h thermal  conductivity,  cal/(cm^) (sec) (°K/cm) 

p density  of  mixture,  g/cnr5 

o collision -diameter,  cm 

cp  equivalence  ratio,  fraction  of  stoichiometric  fuel -oxygen  ratio 

Subscripts : 

eff  effective  value  at  mean  reaction  temperature 
f at  flame  temperature 
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m at  mean  reaction-zone  temperature 

0 initial  conditions 


APPROXIMATE  THEORETICAL  EQUATIONS  USED 
The  Semenov  Bimolecular  Equation 

Zeldovich  and  Frank-Kamenetsky  obtained  an  approximate  solution 
for  the  rate  of  flame  propagation  from  the  differential  equations  for 
heat  conduction  and  fuel  concentration  change  across  the  flame  front- 
This  approximate  solution  was  presented  in  detail  by  Semenov  (refer- 
ence 8) - The  general  approximate  solution  may  be  written  ( see  refer- 
ence 8 , pp-  31,  48,  and  49) : 


, , TTf 

2X  / w dT 


/ n \m  > , m 

_ I ^i\  / X \ 

aopoe:p  (Tf  To^  vV  \Dcpp/f 


(l) 


where,  for  a bimolecular  reaction  between  unlike  molecules 


if  *f 

w dT  » k(a)(b)  dT 

J 0 


RTf 

= aeff  beff  PZ  ' 


-E/RTf 


(2) 


Far  lean  mixtures,  cp<l,  aeff  and  are  computed  by  (fol- 

lowing the  assumptions  of  reference  8,  p-  45)  : 


RT< 


l0  E 


aeff  = “0  T^ 


3eff 


30E 


RTf 


1 - «P\1  - 


Tf ~Toy 


(3) 


For  rich  mixtures,  tp  > 1 : 


aeff  =a0  fT 
f 


RTf 


E 


W 
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beff  ■ 


RTf2 

To_ 
0 Tf  Tf-To 


U) 


For  stoichiometric  mixtures,  equations  (3)  and  (4)  are  identical. 
All  of  the  factors  in  equations  (l)  through  (4)  can,  with  the  exceptions 
of  P and  E,  be  estimated  by  the  extrapolation  of  thermodynamic  tables 
(references  12  and  13)  or  calculated  by  the  relations 
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(Values  calculated  for  combustion-product  mixtures  by  this  simple,  addi- 
tive method  were  within  1 percent  of  values  calculated  by  the  method  of 
reference  14.) 


Cp“^CP,iNi  > 

X = (Cp  + ^ R)  rj/H  (reference  13) 

D «»  1-336  q/p  (reference  15)  \ (5) 


P 

Z 


8xRT 


(M^g) 


fll/2 


M1M2 
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Flame -velocity  predictions  by  the  Semenov  equation  are  evaluated 
herein  as  follows-  Activation  energies  from  low-temperature  reactions 
and  measured  flame  velocities  are  substituted  in  equation  (l)  and  a 
steric  factor  P Is  calculated  for  each  experimental  flame  velocity. 
These  P's  are  then  averaged  to  give  P for  the  group  of  data  points. 
The  ratio  of  the  predicted  flame  velocity  to  the  experimental  velocity 
is  calculated  as  (P/P) l/2 . The  average  deviation  of  the  ratios  from  1 
is  considered  to  be  a measure  of  the  accuracy  of  the  predicted  flame 
velocities.  (The  word  accuracy  is  used  to  denote  the  degree  of  corre- 
spondence between  flame  velocities  calculated  by  the  methods  described 
and  measured  flame  velocities.) 


The  data  for  various  hydrocarbons  consisted  of  maximum  flame 
velocities  (maximum  with  respect  to  equivalence  ratio)  of  56  pure  hydro- 
carbons at  298°  K and  atmospheric  pressure  (references  Z and  3) • These 
data  were  obtained  by  the  NACA  tube  method.  The  hydrocarbons  included 
straight  and  branched  chain  members  of  the  alkane,  alkeae,  alkadiene, 
alkyne,  and  cycloalkane  series  and  bensene-  No  calculations  are  made  by 
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the  Semenov  equation  for  these  data,  because  the  work  required  by  this 
treatment  would  not  be  warranted  in  view  of  the  slight  differences  in 
flame  velocities  observed  for  most  of  the  hydrocarbons  studied.  However, 
it  has  been  shown  by  more  approximate  calculations  (reference  16)  that 
relative  predictions  would  be  satisfactory. 

The  initial -mixture -temperature  data  included  maximum  flame  veloc- 
ities of  methane-air,  propane -air,  and  ethylene-air  mixtures  obtained  by 
a Bunsen  burner  method  (total  area  method,  outside  edge  of  shadow  cast 
by  flame  cone,  reference  5)  at  seven  initial  mixture  temperatures  ranging 
from  200°  to  617°  K.  The  following  activation  energies,  in  kilocalories 
per  gram-mole,  are  reported  in  references  17,  18,  and  19,  respectively: 
methane,  51;  propane,  38;  ethylene,  40.  These  values  are  used  in  the 
Semenov  equation  to  evaluate  the  flame-velocity  predictions-  The  accu- 
racy of  the  predicted  flame  velocities  is  2 to  3 percent.  Specific 
values  sure  tabulated  in  table  I under  "Semenov  Equation". 

Two  kinds  of  mixture-composition  variables  are  studied-  The  first 
is  the  change  of  hydrocarbon  concentration  in  air  over  an  equivalence - 
ratio  range  of  0-7  to  1.3.  Flame-velocity  values  by  the  tube  method  for 
ethylene-air  and  pentane-air  at  298°  K (reference  4) , and  Bunsen  burner 
values  (total  area  method,  outer  edge  of  cone  shadow)  for  methane-air 
at  307°  K and  propane-air  at  302°  K (reference  5)  are  used-  When  the 
effect  of  changes  in  equivalence  ratio  on  the  flame  velocity  is  computed 
by  the  Semenov  equation,  it  is  found  that  the  approximations  used  in 
equations  (3)  and  (4)  are  not  very  consistent  for  the  region  near 
stoichiometric,  particularly  for  equivalence  ratios  between  0-95  and 
1-05-  This  fact  is  shown  by  the  difference  between  the  Semenov  curves 
in  figure  l(a),  where  points  in  the  range  0-95  to  1-05  are  Included  in 
plotting  the  curve,  and  figure  1(b),  where  these  points  are  omitted- 
If  this  critical  region  is  emitted,  the  accuracy  of  the  predicted  flame 
velocities  over  the  equivalence  range  is  3 to  7 percent-  Specific 
deviations  are  listed  in  table  I- 

The  second  composition  variable  studied  is  the  ratio  of  oxygen  to 
oxygen  plus  nitrogen  a * OgAOg+Ng)  • Maximum  flame  velocities  for 

several  a ratios  were  available-  for  the  following  fuels  and  ranges  of 
a:  propane,  0.17  to  0.50;  ethylene,  0.17  to  0-35;  and 

2,2,4-trimethylpentane,  0.21  to  0-50.  The  accuracy  of  prediction  of 
these  data  is  3 to  17  percent  (table  I) . 

It  may  be  noted  in  table  I that  flame-velocity  predictions  by  the 
Semenov  equation  for  ethylene  flames  are  accurate  to  4 percent  for 
variations  in  both  composition  and  initial  tenperature . For  these 
predictions  (as  discussed  previously),  each  type  of  data  was  considered 
separately,  and  steric  factors  were  calculated  from  the  data.  Now 
l usually,  when  it  is  desired  to  predict  flame  velocities,  only  one  type 
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of  data  is  available  - for  example,  fuel-air  ratio  data  at  atmospheric 
pressure  and  room  temperature.  It  is  therefore  Important  to  know  how 
well  the  effects  of  initial  temperature  and  oxygen  concentration  can  be 
predicted  from  such  data.  The  ethylene  data  are  used  as  an  example- 

The  steric  factor  P calculated  from  the  equivalence -ratio  data 
and  the  low -temperature  activation  energy  of  40  kilocalories  per  gram- 
mole  are  used  in  the  Semenov  equation  to  predict  flame  velocities  over 
the  ranges  of  initial  temperature  and  oxygen  concentration  covered  by 
the  experimental  data.  In  both  cases  the  predicted  flame  velocities 
deviated  from  the  normalized  measured  flame  velocities  by  an  average  of 
+14  percent,  the  maximum  deviation  being  +24  percent  for  the  oxygen- 
concentration  data.  (Because  different  methods  of  measuring  flame 
velocity  were  used,  differences  in  values  for  mixtures  of  the  same  com- 
position and  temperature  were  found.  For  these  calculations  all  flame 
velocities  were  normalized  by  a simple  ratio  factor  to  the  values  for 
the  initial- temperature  data.) 


Tanford-Pease  Equation 

Tanford  and  Pease  (references  9 and  20)  equated  the  amount  of 
product  formed  in  the  combustion  zone  by  a second-order  reaction  between 
fuel  molecules  and  hydrogen  atoms  (or  other  active  particles  such  as 
hydroxyl  radicals  or  oxygen  atoms)  to  the  amount  of  product  formed  at 
the  flame  front  by  conversion  of  the  fresh  gas  expressed  in  terms  of 
initial  conditions  and  flame  velocity-  An  approximate  solution  for  flame 
velocity  from  this  equation  was  given: 


Three  methods  of  evaluating  the  predictions  of  the  Tanford  and 
Pease  equation  are  used  herein-  For  all  three,  the  following  calcula- 
tions are  the  same:  (a)  flame  temperature  Tf  and  active  particle 

concentrations  -p^  are  calculated  assuming  adiabatic  thermal  equilibrium 
by  a matrix  method  (reference  21);  (b)  the  mean  combustion-zone  tempera- 
ture is  assumed  to  be  0.7  Tf  (reference  9);  (c)  diffusion  coefficients 
D^,  are  calculated  from  Dm  - Df  (0-7  Tf/To)^‘°7,  where  Di  is  the 
diffusion  coefficient  at  intial  temperature  calculated  by  the  Stefan- 
Maxwell  equation  (reference  22);  (d)  the  reconfloination  factor  Bj  is 
calculated  by  the  method  of  Tanford  (reference  10)  for  the  hydrogen 
atoms  and  is  assumed  to  be  one  for  OH  and  0;  (e)  the  ratio  LmQ'n/Q  is 

calculated  from  a knowledge  of  the  over-all  oxidation  process  and  the 
initial  concentrations  of  reactants. 


...  •/ 
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The  three  methods  of  evaluation  differ  in  the  calculation  of  k*. 
For  the  first  method  only  one  chain  carrier,  the  hydrogen  atom,  is  con- 
sidered- For  this  case  % values  are  calculated  from  single-point 
flame-velocity  determinations  by  equation  (6).  For  the  second  method, 
II,  0,  and  OH  are  considered  to  be  the  chain  carriers  and 
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Single-point  ka  values  are  calculated.  For  these  two  methods,  the 
ratio  of  the  predicted  flame  velocity  to  the  experimental  flame  velocity 
is  calculated  as  (lt^/k^)1/2  and  the  average  deviation  of  the  ratio 

from  1 is  used  as  a measure  of  the  accuracy  of  the  predictions. 

For  the  third  method,  an  Arrhenius  type  temperature  dependence  of 
the  rate  constant  is  used: 

- PaZa  exp  (-EjK(0.7  Tf))  (8) 

The  following  method  is  used  to  calculate  an  activation  energy  Ea  for 
best  fit  from  the  initial-temperature  data.  Several  values  of  £a  are 
substituted  in  the  equation  and  for  each  a value  of  (P/F)1/^  is  cal- 
culated. Then  6,  the  percentage  mean  deviation  in  (P/PW2  fo£  the 
group,  is  plotted  against  the  Ea  used.  The  "fitted”  value  of  Ea  is 
the  one  which  gives  the  minimum  &•  All  other  calculations  by  the  third 
method  are  made  by  equations  (7)  and  (8),  measured  flame  velocities,  and 
Ea  from  initial -temperature  data.  The  accuracy  is  calculated  as  the 
mean  deviation  of  (P/P)1?2  from  1- 

The  predictions  of  the  Tanfard-Pease  equation  are  evaluated  by 
methods  1 and  2 for  the  maximum  flame  velocities  of  a series  of  hydro- 
carbons (reference  23) - Table  I shows  that  the  accuracy  of  prediction 
of  flame  velocity  is  3-3  percent  considering  hydrogen  atoms  (kg)  only 
and  3-5  percent  considering  all  three  active  particles  (1%).  These 
mean  deviations  are  small  enough  to  indicate  either  that  the  "rate  con- 
stants" are  not  temperature-dependent,  in  which  case  the  steric  factors 

are  of  the  order  of  10'3,  or  that 

stant  for  the  hydrocarbons  studied,  which  is  probable.  If  the  activation 
energy  is  of  the  order  of  10  kilocalories  per  gram-mole,  the  steric  fac- 
tor is  of  the  order  of  10~^ . 

All  three  methods  of  evaluation  are  used  with  the  initial- 
temperature  data-  The  accuracy  of  the  flame  velocity  predictions  is 


|exp  (r ( 0 - 7 Tf)) 


T^  is  nearly  con- 
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9 to  12  percent  (table  I)  for  the  first  method  (kg)  and  6 to  8 percent 
for  the  second  method.  The  use  of  method  3,  which  includes  the  tempera- 
ture dependence  of  the  rate  constant,  improves  the  prediction  so  that 
the  accuracy  is  about  2 percent.  The  calculated  activation  energies  and 
steric  factors  are:  for  methane,  17  and  0-19;  and  for  ethylene,  18  and 

0.21*  The  calculated  activation  energy  for  ethylene  is  high  when  com- 
pared with  the  experimental  activation  energy  for  the  reaction 
H + C2H4— >C2Hs,  which  is  2-6  kilocalories  (calculated  from  data  in 
reference  24) - The  relative  predictions^  axe  not  appreciably  changed  by 
assuming  Bg  ■=  1 or  by  assailing  Dmrc  Tn  which  was  used  in  reference  7 
instead  of  E^coq/p. 

If  the  entire  equivalence -ratio  ranges  which  were  covered  experi- 
mentally by  the  flame -velocity  data  are  used  for  the  prediction  of 
velocities  by  the  Tanford-Peuse  method,  the  accuracy  of  prediction  is 
very  poor  as  illustrated  in  figure  1(a)  for  ethylene-  If  the  equivalence 
range  is  limited  to  one  near  the  maximum  flame  velocity  (for  example, 
cp  = the  accuracy  of  the  predictions  is  I to  6 percent 

(table  I)  when  all  active  particles  are  considered.  However,  if  only 
hydrogen  atoms  are  used  the  accuracy  is  very  poor  - 17  to  65  percent. 

The  variations  of  flame  velocity  with  oxygen  concentration  were 
calculated  for  all  active  particles.  The  accuracy  of  prediction  is  4 to 
6 percent. 

The  ethylene  data  are_again  used  as  an  example  to  show  how  well  an 
average  "rate  constant”  (1%  from  method  2)  from  the  fuel-air  ratio 
data  predicts  the  effects  of  temperature  and  oxygen  concentration  on  the 
maximum  flame  velocity.  Predicted  flame  velocities  for  various  tempera- 
tures differ  from  the  normalized  experimental  values  by  -9  to  +13  percent 
with  an  average  deviation  of  6 percent.  Predicted  flame  velocities  for 
various  oxygen  concentrations  differ  from  experimental  values  by  +8  to 
+22  percent  with  an  average  deviation  of  12  percent. 


The  'Manson  .Equation 


The  momentum  relationship  between  the  flame  (or  combustion  wave) 
velocity  and  the  pressure  drop  across  a plane,  steady-state  flame  front 
may  be  expressed  in  the  form  (reference  25,  pp-  241): 


U » 


pQ  P0-Pf, 


(9) 
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Manson  (reference  11)  suggested  that  the  small  pressure  drop  could  he 
caused  by  the  projection  of  hydrogen  atoms  into  the  unburned  gas.  Because 
the  H atoms  would  recombine  to  Hg  at  the  unburned  gas  temperature, 

Ap  was  assumed  to  be  one-half  the  equilibrium  H atom  pressure  reduced 
to  unburned  gas  temperature-  In  the  present  paper,  Ap  is  calculated  as 


Tr 


both 


PH 


and 


doh 

P0H5J- 


DoNTo 

Po  »hK' 


The  second  method  of 


calculating  Ap  takes  into  account  three  of  the  lightest  and  most 
abundant  species  which  might  be  considered  active  particles  and  assumes 
that  all  of  these  would  recombine  in  pairs.  The  average  proportionality 
factor  for  a group  of  data  is  calculated  by  dividing  the  experimental 
flame  velocity  by  the  right-hand  side  of  equation  (9)  for  each  point  and 
averaging  the  quotients-  These  average  empirical  proportionality  factors 
are  designated  % for  the  first  method  and  Kj  for  the  second-  The 


accuracy  of  the  predicted  flame  velocities  is  considered  to  be  the 
average  deviation  of  %/%  or  %/Kj  from  1-  For  the  maximum  flame 


velocities  of  the  different  hydrocarbons,  the  accuracy  of  prediction  by 
the  Manson  equation  for  either  method  of  evaluation  is  about  4 percent 
(table  I).  For  the  first  method  (Ky),  the  proportionality  factor  is 
0-5;  that  is,  the  predicted  value  is.  twice  the  experimental  value- 


For  the  initial-temperature  data  the  predicted  flame  velocities 
differ  from  the  experimental  velocities  by  10  to  14  percent.  When  the 
total  range  of  the  equivalence -ratio  data  is  considered,  the  accuracy  of 
prediction  of  flame  velocities  by  the  Manson  equation  and  either  of  the 
methods  of  evaluating  Ap  is  poor-  If  the  equivalence  range  is  limited 
as  for  the  Tanford-Pease  method  and  with  all  active  particles,  predic- 
tions accurate  to  2 to  4 percent  are  obtained  (table  I) - ■ The  accuracy 
of  prediction  of  flame  velocities  for  mixtures  containing  different 
concentrations  of  oxygen  varies  from  hydrocarbon  to  hydrocarbon.  For  - 
2,2,4-trlmethylpentane,  the  accuracy  is  about  1 percent;  for  ethylene, 

5 percent;  and  for  propane,  15  percent. 


The  ethylene  data  are  again  considered  in  order  to  determine  the 
accuracy  with  which  the'Kj-  from  the  equivalence -ratio  data  predicts 
the  effects  of  temperature  and  oxygen  concentration  on  flame  velocity. 

The  deviations  for  the  initial-temperature  data  varied  from  -3  to  +45  per- 
cent with  a mean  deviation  of  27  percent.  For  the  oxygen  data  the 
deviations  varied  from  +1  to  +19  percent  with  a mean  value  of  +10  percent- 


Empirical  Correlations 

A correlation  was  obtained  (reference  26)  by  which  maximum  flame 
velocities  of  various  hydrocarbons  in  air  could  be  predicted  with  an 
average  deviation  of  2 percent.  The  predicted  flame  velocity  is  calcu- 
lated from  the  sum  of  the  contributions  of  various  H-C  bonds  in  the  fuel 
molecule  according  to  the  following  relationship: 
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U - NaKa  + %%  + NcKc  + • • • 

where  NC,  UD,  NE,  Np,  N0,  and  Ng  are  the  numbers  of  methane, 

primary,  secondary,  tertiary,  alkene,  alkyne,  cyclohexyl  and  aromatic 
C-H  bonds,  respectively,  per  unit  volume  of  hydrocarbon -air  mixture, 
and  Ka,  Kg,  Kg,  . . . are  the  flame-speed  coefficients  of  these  bonds. 

For  the  special  cases  of  C-H  bonds  on  carbon  atoms  placed  alpha  to  the 
alicyne  C = C bond,  a factor  M = 0.96  was  introduced  into  terms  repre- 
senting these  alpha  bonds-  Correlation  coefficients  (table  II)  estab- 
lished from  34  hydrocarbons  excluding  ethylene  then  gave  an  average 
deviation  of  1.9  percent  in  the  ratio  of  predicted  to  measured  flame 
velocity. 

A correlation  was  also  observed  between  the  equivalence  ratio  for 
maximum  flame  velocity  and  the  total  bond  dissociation  energy  of  the 
fuel-  As  shown  in  table  III,  the  total  bond  dissociation  energy  per 
unit  volume  of  the  hydrocarbon-air  mixture  corresponding  to  the  maximum 
flame  velocity  Is  nearly  constant,  with  an  average  deviation  of  only 
0.9  percent  from  the  average  value.  This  deviation  is  less  than  one- 
third  of  that  obtained  by  simply  assuming  that  the  maximum  will  occur 
at  an  average  equivalence  ratio  of  1.15.  Bond  energy  calculations  for 
37  hydrocarbons  similar  to  those  presented  in  table  III  had  an  average 
deviation  of  1.6  percent- 

For  the  equivalence -ratio  data,  linear  relationships  were  found 
between  flame  velocity  and  the  logarithm  of  the  equivalence  ratio  cp, 
where  cp  <cpmnY-  An  interesting  variation  is  the  plot  of  flame  velocity 

against  the  total  bond  dissociation  energy  of  the  fuel  per  unit  volume 
of  mixture  in  figure  2 for  pentane,  ethylene,  and  propyne  mixtures  with 
air.  It  may  be  seen  that  the  correlation  is  linear  for  the  lean  mix- 
tures and  that  the  line  extrapolates  to  the  lean  limit  for  flame  propa- 
gation (where  U = 0)  as  determined  in  the  1-inch  flame  tube-  For  these 
data,  determined  by  the  tube  method  for  three  hydrocarbons , good  esti- 
mates of  the  maximum  flame  velocities  are  obtained  by  reading  from  the 
curves  at  a total  bond  dissociation  energy  of  1-75  kilocalories  per 
liter.  This  constant  bond-energy  value  of  1-75  kilocalories  per  liter  for 
this  correlation  does  not  equal  the  "bond  energy  of  the  hydrocarbon  con- 
centration corresponding  to  the  maximum  flame  velocity  which  was  pre- 
viously calculated  to  be  1-96  kilocalories  per  liter.  Work  with  other 
flame-velocity  data  shows  that  both  the  position  of  the  horizontal  line 
and  the  slopes  of  the  lineB  for  specific  hydrocarbons  depend  somewhat 
on  the  method  of  measurement  of  flame  velocity. 

Empirical  equations  predicting  the  effect  of  initial  mixture  tem- 
perature on  flame  velocity  with  mean  deviations  of  approximately  2 per- 
cent for  the  range  from  200°  to  615°  K are  (reference  5) : 
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for  methane 

-5  2.23 

U - 10  + 7.40  X 10  ° T0 

for  propane 

U « 10  + 3-42  X10"4  T02'°° 

for  ethylene 

U =»  10  + 2.59  X10'3  Tq1'74 

Empirical  relationships  were  found  which  predict  the  effect  of 
oxygen  concentration  (for  the  experimental  range  covered)  and  limited 
temperature  (311°-422°  K)  on  the  flame  velocity  for  the  three  hydro- 
carbons studied  (reference  6) : 

for  2,2, 4-trimethylpentane 

U = 0-133  T01‘4°(a-0.120) 


for  propane 

U = 0-766  T0I,i6(a-0.133) 


for  ethylene 

U =■  0-998  T01'18(a-0.133) 

The  accuracy  of  prediction  was  about  3 percent  for  ethylene  and 
2, 2, 4-trimethylpentane  and  6 percent  for  propane- 


CONCLUDING  REMARKS 

The  semitheoretical  calculations  of  this  paper  are  dependent  on  the 
assumption  that  thermal  equilibrium  is  attained  at  the  end  of  the  com- 
bustion zone.  Actual  flame  temperatures  and  product  concentrations 
would  differ  from  equilibrium  values  if  (a)  equilibrium  is  not  attained, 

(b)  any  appreciable  reaction  takes  place  at  the  initial  temperature,  or 

(c)  there  is  a large  chain  branching  term  affecting  the  radical  distri- 
bution in  and  ahead  of  the  flame  zone-  These  calculations  are  also 
dependent  on  the  Inherent  assumptions  of  the  approximate  theoretical 
equations  and  the  methods  used  to  calculate  transport  properties-  Within 
these  limitations,  the  results  of  the  calculations  may  be  generalized 

as  follows : 
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1-  Tie  Uemenov  bimolecular  equation  (thermal  mechanism),  the 
Tanford-P jasu  square -root  law  (active -particle  diffusion  mechanism),  and 
the  Manse  a modification  of  the  momentum-pressure -drop  equation  (including 
diffusion  of  active  particles,  hut  not  chemical  kinetics)  will  all  predict 
relative  changes  in  flame  velocity  caused  by  changes  in  hydrocarbon 
structure,  initial  temperature,  equivalence  ratio  ( limited  range) , or 
oxygen  cc nces.tration.  The  accuracies  of  the  predictions  are  summarized 
in  table  I. - * 

2.  'vOw-1  emperature  activation  energies  may  be  used  in  the  Semenov 
equation  to  give  good  relative  predictions  of  flame  velocities  for  var- 
iations '.\i  initial  temperature,  equivalence  ratio,  and  oxygen  concentra- 
tion. Only  the  bimolecular  fuel-oxygen  reaction  is  considered  in  this 
paper;  ether  lombinations  such  as  fuel  and  hydrogen  atoms  should  be 
investigated  to  determine  whether  the  kinetics  indicated  for  such  other 
reaction..;  are  not  more  probable  than  for  the  fuel-oxygen  case. 

3 .he  average  "specific  rate  constants"  (k^)  obtained  with  the 
Tanf or  i -Pease  equation  for  molecular  fuel  - active  particle  reactions 
behave  ts  Arruenius  rate  constants  only  with  the  data  on  the  effect  of 
initial  temperature,  not  with  the  equivalence -ratio  or  oxygen- concentra- 
tion k* ta.  The  relative  predictions  obtained  by  this  equation  are  not 
very  nsitiv--  either  to  the  temperature  dependence  assigned  to  the 
diffusion  coefficients  or  to  the  recombination  factor  computed  for  hydro- 
gen aura;;  thi  ne  factors  do  appreciably  affect  absolute  predictions  of 
flame  velocity  or,  conversely,  determinations  of  steric  factors  from 
exper L: usntal  l lame  velocities-  Regarding  active  particles  considered, 
better  rt cults  are  generally  obtained  for  hydrocarbon  flames  when  H, 

OH,  e:ii  0 are  considered  rather  than  H alone;  this  is  particularly  true 
for  studies  of  the  effect  of  equivalence  ratio. 

4.  For  the  equivalence -ratio  and  oxygen-concentration  data,  con- 
siderably better  relative  predictions  are  obtained  from  the  Manson  equa- 
tion v'hen  the  pressure  drop  across  the  flame  front  is  considered  to  be 
due  U.  H , UH,  and  0 than' when  only  H atoms  are  taken  into  account- 
Furt  er  thought  should  be  given  as  to  how  this  pressure  drop  should 
be  craputed,  taking  into  account  any  net  diffusional  flow  of  any 
comprnent  between  the  flame  front  and  the  unbumed  gas. 

3.  The  example  presented  by  the  ethylene  data  shows  that  an  empir- 
ically determined  rate  constant  or  proportionality  factor  from 

equl v n le no e -ratio  data  at  atmospheric  pressure  and  room  temperature  may 
pred'et  the  effect  of  initial  temperature  or % oxygen  concentration  within 
approximately  20  percent  by  the  Semenov  or  Tanford-Pease  equation  and 
within  approximately  30  percent  by  the  Manson  equation. 

6-  For  engineering  applications,  the  effects  of  the  parameters 
s tut.  led  could  be  estimated  Just  as  satisfactorily,  and  more  easily,  by 
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TABLE  I - PERCENTAGE  MEAN  DEVIATIONS  IN  RATIO  OF  PREDICTED  TO  MEASURED  FLAME  VELOCITY 
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at 


t/ta  determined  by  tube  method  at  room  temperature  (reference  2) . 

tta  determined  by  Bunsen  burner  method,  shadow  cone  (reference  5) • 

ita  determined  by  Bunsen  burner  method,  schlieren  cone  (references  6 and  7). 
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Single-bond  energies  from  reference  27.  Multiple -bond  energies  from 
reference  28:  151-2  kcal/g-mole  for  C=C  in  olefins  and  cyclic  compounds 

198-5  kcal/g-mole  for  CSC  bond- 


TABLE  XI  - EMPIRICAL  COEFFICIENTS  FOR  CALCULATING  MAXIMUM  FLAME 
VELOCm  FROM  HYDROCARBON  STRUCTURE  (REFERENCE  26) 


Type  C-H  bond 

Coefficient 

Methane 

ka 

35-2X10'19 

Primary 

kb 

42.5 

Secondary 

«C 

47.5 

Tertiary 

Kd 

45.4 

Cyclohexyl 

kg 

50-5 

Alkene 

% 

80.7 

Aromatic 

% 

84-3 

Alkyne 

KF 

223-9 

TABLE  III  - BOND  DISSOCIATION  ENERGIES  FOR  C4  AND  Cg  HYDROCARBONS 
AT  CONCENTRATION  FOR  MAXIMUM  FLAME  VELOCITY 


Hydrocarbon 

Equiva- 

lence 

ratio 

[Deviation  Dissocia- 

°f  tPmay  tion  energy 

from  average  of  hydro- 

Dissocia- 
tion energy 
of  mixture 

Devia- 

tion 

from 

j 

for  max- 

j^iax 

carbon 

at 

average 

lmum  U, 

/ kcal  \ 

1 kcal  \ 

(percent) 

S (percent) 

L— .. 

' u>  I 

\ XI  / 

Butane 

1-09 

5-5 

1250 

1-93 

2.0 

Hexane 

1.16 

0.9 

1791 

2.00 

1.5 

2 -Mettaylpropane 

1-11 

3-6 

1257 

1-94 

1.5 

Z, 2-Dimethyl- 

butane 

1.12 

2-5 

1821 

1-97  - 

0.0 

Butene-1 

1-17 

1.7 

1133 

1.98 

0.5 

Hexene -1 

1.16 

0.9 

1680 

2-02 

2.5 

Isobutene 

1-14 

0-9 

1147 

1-97 

0.0 

2 -Methylpentene - 1 j 

1-19 

3-4 

1603 

1-99 

1-0 

Butyne-1 

1-17 

1.7 

1016 

1-95 

1.0 

Hexyne-1 

1-21 

5.0 

1564 

1-97 

0-0 

4 -Methylpentyne 

1.18 

2-5 

1574 

1-95 

1.0 

Benzene 

1.34 

14.2 

1305 

1.96 

0.5 

Average 

1.15 

3-29 

— 

1-96 

0.9 

Burning  velocity,  ca/sec 


Figure  2.  - Empirical  correlation  for  equivalence  ratio  data 


